Coherent control of nanoscale light localization in metamaterial: creating 
and positioning a sub-wavelength energy hot-spot 
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Precise control and manipulation of optical fields on a nanoscale is one of the most important and 
challenging problems in "nanophotonics" . Since optical wavelength is on a much larger microscale, it 
is impossible to employ conventional focusing for that purpose. We show the strong optically-induced 
interactions between discrete meta-molecules in a metamaterial system and coherent monochromatic 
continuous light beam with a spatially-tailored phase profile can be used to prepare a sub- wavelength 
scale energy localization. Well isolated energy hot-spots as small as A/10 can be created and 
positioned at will on the metamaterial landscape offering new opportunities for data storage and 
imaging applications. 



PACS numbers: 

In 2002 a method was suggested, which was 
based on tailoring phase modulation of ultrashort 
optical pulses in the time domain to achieve coher- 
ent control of spatial distribution of the excitation 
energy in complex inhomogeneous nanosystems 
[l|. In such systems the localized plasmon fre- 
quencies vary from nanoscale feature to nanoscale 
feature and thus can be correlated with their po- 




FIG. 1: Coherent control of light localization in ar- 
ray of interacting meta-molecules: artistic impression. 
A light beam with modulated phase profile excites a 
metamaterial array of interacting plasmonic resonators 
leading to the formation of the energy hot-spot at the 
metamaterial. This spot may be moved from a meta- 
material cell to another metamaterial cell by tailoring 
the phase profile of the incident beam. 



sitions. The pulse phase modulation will cause 
the exciting field to take energy away from surface 
plasmons localized in those parts of the system 
where the oscillations are out of phase with the 
driving pulse and move it, with time, to the sur- 
face plasmon excitations in other parts where such 
oscillations occur in phase with the driving pulse. 
Later a similar idea based on tailoring of ultra- 
fast pulses in the vicinity of silver nanostructures 
through adaptive polarization shaping has been 
successfully demonstrated @, Illuminating a 
sub-wavelength diffraction grating by an ampli- 
tude and phase modulated continuous-wave light 
beam has been suggested as a method to generate 
sub-diffraction light spots 

Here we show that strong interactions between 
regularly spaced plasmonic resonators in a pla- 
nar metamaterial, a two-dimensional array of 
nanoscale meta-molecules, can result in a sub- 
wavelength localization of optical energy. A de- 
sired nanoscale light hot-spot can be engineered 
simply by adjusting the far-field spatial profile of 
an incoming monochromatic coherent continuous 
light beam. We show by using a simple model 
how the localization results from a co-operative re- 
sponse of the metamaterial array of altimetrically 
split ring plasmonic resonators, providing physi- 
cal insight into the energy localization effect. The 
essential requirement in the sub- wavelength local- 
ization process of the incoming wave is that in- 
dividual meta-molecules act as closely-spaced dis- 
crete scatterers whose strong inter-meta-molecular 
interactions modify the energy eigenspectrum of 
the system. This is in contrast with the homog- 
cnization approach to electromagnetic properties 
of metamaterials treated as a continuous medium 
instead of a collection of discrete meta-molecules. 
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We show that evanescent field energy can be 
localized almost solely at a single meta-molecule 
while a group of neighboring meta-molecules re- 
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main unexcited. Now by simply adjusting the 
spatial phase profile of the incident beam the 
nanoscale hot-spot can be moved at will from one 
meta-molecule to another, providing an efficient 
technique for a sub-wavelength scale optical con- 
trol and manipulation of a metamaterial system. 
In contrast with localization techniques relying on 
laser-pulse excitation [IHl] our method is based on 
a continuous wave source and does not depend on 
a transient redistribution of energy between nano- 
objects. This makes the localization in a metama- 
terial system much simpler to implement. More- 
over, it does not require the nanoscale system to 
be spatially inhomogeneous (such as a rough sur- 
face) and works with periodic, regular planar array 
of identical nano-objccts opening the opportuni- 
ties for microscopy and data storage applications 
which will be discussed below. 

The idea of coupled resonators illuminated by 
a phase-modulated beam has a clear mechanical 
analogy in two coupled identical oscillators (mod- 
eling nano-objects) that may be driven to vastly 
different amplitudes by setting up a phase delay 
between otherwise identical coherent mechanical 
forces driving them. A system consisting of reso- 
nant wires has previously been used as a metalens, 
transferring sub-wavelength features of a evanes- 
cent field to the spectral properties of propagat- 
ing microwaves at different frequencies that are 
detectable in the far-field [5[ (see also Ref. Q). 

We demonstrate that controllable nanoscale lo- 
calization of optical energy can be achieved at 
an artificial nanostructured metamaterial surface, 
consisting of an array of altimetrically split ring 
plasmonic resonators, when such structure is il- 
luminated with a coherent continuous light wave 
(a laser) those spatial variation of the phase at 
the metamaterial plane can be tailored. The re- 
quired spatial variation of phase are slow and shall 
take place on the scale longer than the wavelength 
of light. This can be routinely achieved experi- 
mentally using widely commercially available liq- 
uid crystal spatial light modulators. The choice 
of asymmetrically split ring array as the nanoscale 
landscape was motivated by the recent studies that 
established a key role of inter-meta-molecular in- 
teraction in forming the electromagnetic response 
of such metamaterials @, H|. We will see below 
that the existence of strong interactions between 
plasmonic resonators is the key requirement for 
nanoscale localization. Figure [1] shows artistic rep- 
resentation of the metamaterial array excited by 
the wave with a tailored phase profile creating a 
hot-spot at the other side of the array, at its im- 
mediate proximity. 

Our analysis of this system is based on di- 
rect quantitative numerical simulation using a 3D 
Maxwell solver illustrated by qualitative model- 
ing the metamaterial as a regular array of dis- 
crete interacting dipolar scatterers. In all numeri- 
cal simulations we concentrated on the visible and 



the near-field parts of the spectrum, assumed that 
metamaterial structures are constructed of gold 
plasmonic resonators and used realistic material 
parameters and Joule loss factors by using the well 
established data for the dielectric parameters of 
the metal Q. 

We demonstrate the idea of coherent control 
of light localization using the excitation beam 
with a periodic phase modulation across its wave- 
front: even this most simple profile leads to a 
profound light localization at the array of inter- 
acting meta-molecules at certain frequencies of 
the driving field. Fig. [5] shows intensity maps 
at the metamaterial array driven by a monochro- 
matic field manipulated such that, upon arrival 
at the metamaterial surface, the electric field has 
constant amplitude, a fixed polarization along y- 
direction, and a spatially varying phase 4>{x, y) = 
(A<fi max /2) sin(«;x) sin(Ky), where k — 2ir/(6a) 
and a is the spacing between adjacent meta- 
molcculcs. This corresponds to the incident wave- 
front modulated along x and y direction with a 
period of six unit cells of the metamaterial array. 
The array is a gold film of 50 nm thickness per- 
forated with 25 nm slits of split square shape, as 
presented in Fig.^h). The metamaterial unit cell 
is 200 nm x 200 nm. The metamaterial nanostruc- 
ture has a well defined plasmonic absorption res- 
onance at 880 nm corresponding to a local maxi- 
mum in transmission and a minimum in reflection 
[Toj . We investigated field localization in the im- 
mediate proximity to the array at this wavelength 
of excitation and at two other wavelengths above 
(1150 nm) and below (700 nm) the resonance, at 
local minima of absorption and transmission. 

The main characteristic features of the coherent 
control process are illustrated on Fig. [2] Fig.^b-c- 
d-e) show the intensity distribution at the array as 
a function of the amplitudes of spatial phase mod- 
ulation. As expected, for a plane wave front of the 
incident beam (A<fi max = 0) all meta-molecules in 
the array are excited equally. With an increase of 
the modulation depth, energy tends to concentrate 
at two opposite quadrants of the 6x6 section of the 
array. At (A(f> rnax — it) two hot-spots are formed 
which are localized at a single meta-molecule six 
unit cells apart from a neighboring hot-spot along 
both x and y directions. For a fixed amplitude 
of phase modulation the localization is most pro- 
nounced at the plasmonic absorption resonance at 
880 nm: de-tuning away from the resonance de- 
stroys the localization, as illustrated on Fig. [21 e- 
f-g). Fig. EJa) shows intensity cross-section of 
the resonant hot-spot (880 nm) which has a foot- 
print (at half maximum) of only 60 nm x 133 nm 
(0.0078 ^m 2 ). 

The main features of the coherent control may 
be illustrated by a simple model in which we treat 
current oscillations in individual meta-molecules 
as point-like dipoles. An incident driving field 
of a given wavelength excites the current oscilla- 
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FIG. 2: Coherent control of energy localization at metamaterial array of interacting plasmonic resonators. Colour 
maps (a)-(d) show the build-up of light localization with the increase of sinusoidal phase modulation of the driving 
field, A(f>. Note two intensity peaks formed at the opposite quadrants of the landscape at resonance conditions at 
A 880 nm and A0 = n. Maps (d)-(f) show the light localization effect at different wavelength for A(/> = it: tuning 
the excitation wavelength away from the resonance (e) A = 700 nm and (f) A = 1150 nm destroys the localization. 
Figure (g) show the far-field transmission, reflection and absorbtion spectra of the metamaterial array and its 
unit cell (h). Localization is most pronounced at A = 880 nm corresponding to the absorbtion resonance 77. 



tions, producing scattered radiation. These emit- 
ted fields in turn drive the corresponding oscilla- 
tions in neighboring meta-molecules. The system 
dynamics can be described by the monochromatic 
version of Maxwell's wave equations, for instance, 

(V 2 + k 2 )D+ = -V x (V x P+) - ikV x M+ (1) 

where D + (f,t), P + (f,t) and M + (f,t) are the 
positive frequency components of the electric dis- 
placement, polarization, and magnetization den- 
sities, respectively. For an array meta-molecules 
labelled by an index j, at position Tj, with elec- 
tric dipole moments dj (t) and magnetic dipole mo- 
ments jUj(t), the polarization and magnetization 
densities are given by P(r,t) — J^V dj(t)5(r — fj) 

and M(f,t) — J2j rj), respectively. This 

model accounts for short range quasi-static dipole- 
dipole interactions between meta-molecules and 



long range interactions mediated by the radiated 
field and radiation losses. As a result, the mate- 
rial displays collective modes of oscillation, each 
with a distinct resonance frequency and collective 
radiative damping rate. The majority of collec- 
tive modes do not couple directly to the uniform 
driving field because only collective modes with 
translational symmetry can couple to such a field. 
However, by allowing spatial variation in the driv- 
ing field (for instance a spatial phase modulation, 
as illustrated above), one shall be able to excite a 
superposition of collective modes. A superposition 
of such modal excitations can lead to the forma- 
tion of hot-spots located only at a few isolated 
meta-molecules. 

The effect of strong inter-meta-molecular inter- 
actions in the energy localization may be illus- 
trated by a simple example of two electric dipoles 
with a separation much less than the wavelength 



4 









X FWHM~ 133nm 

(-0.15 A) j 

Y FWHM~ 60nm I 

(-0.07 A) 










/ 


: 



440 880 1320 

X/Y distance, nm 




FIG. 3: The main cross-sections of the energy hot- 
spot (a) created at the metamaterial landscape (see 
Fig. 2d) indicate a sub- wavelength localization with the 
total footprint at the peak of about 1 percent of A 2 . 
Colour maps (b) and (c) show the intensity and phase 
of excitations within the model of individual interact- 
ing dipoles representing the met a- molecules of array. 

of light. The collective state of the dipoles ex- 
hibits super-radiant (enhanced radiative damping 
rate) and sub-radiant (suppressed radiative damp- 
ing rate) eigenmodes. In the first eigenmode the 
two dipoles are in-phase and in the second one they 
are it out-of-phase. If we could prepare an incident 
field that excites an equal superposition of the two 
modes, it could then localize the energy in one of 
the dipoles. The crucial part in the emergence 
of the sub-radiant and super-radiant states in a 
two-dipole system are recurrent scattering events 
[Til Il2l | - processes in which a wave is scattered 
more than once by the same dipole - which cannot 
be described by the standard continuous medium 
electrodynamics, necessitating a discrete scatterer 
model of the metamaterial system. We present a 
less trivial example in the supplementary material. 
A numerical example of field localization achieved 
within the dipole model and showing all charac- 
teristic features of the coherent process control is 
illustrated in Fig. Ob). 

The opportunity to localize energy hot-spot at 
a single isolated meta-molecule provides for in- 



teresting applications. Indeed, by simply shift- 
ing phase profile function of the driving beam 
4>(x, y) = (A(f> max /2) sm(nx + Sx) sin(/-q/ + Sy) on 
5x and Sy, integers of the lattice period, one can 
move the hot-spot at will across the metamaterial 
landscape. In what has been illustrated above, the 
hot-spots can be spaced by at least six unit cells 
from one another and thus moved at will in this 
space. Our computational capabilities do not al- 
low for numerical modeling of driving fields mod- 
ulated at a larger scale which would involve mod- 
eling of a larger array of meta-molecule. We ar- 
gue, however, that similar localization effects could 
take place in larger arrays where individual hot- 
spots are localized in much larger dim areas. More- 
over, we argue that localization will be improved 
by appropriate selection of metamaterial and by 
controlling losses (for instance by using silver in- 
stead of gold) that will ensure long-range interac- 
tions between meta-molecules. 

Controlling the hot-spot position on the meta- 
material chess-board may be used for high-density 
data storage application in particular with phase- 
change nano-particles 13] placed at the meta- 
molecules and switched at the hot-spot. Moreover, 
the light localization technique can be used in a 
new type of optical imaging instruments with sub- 
wavelength resolution: one can envisage that an 
object of study (say, a cell) is placed on the meta- 
material surface with a light-collecting optics and 
a photo detector above it. Moving the hot-spot 
from one place to another shall lead to the hot-spot 
being exposed to different scattering/absorption 
regimes which will affect intensity at the photo 
detector. A map of such intensities can be build 
to reconstruct a sub-wavelength image of the ob- 
ject. Here the metamaterials lattice pitch defines 
imaging resolution. 

In conclusion, we demonstrated coherent con- 
trol of light localization in photonic metamaterial. 
Well isolated energy hot-spots as small as A/10 can 
be created and positioned at will on the plasmonic 
landscape by illuminating an array of interacting 
plasmonic resonators with a coherent continuous 
light beam with a tailored spatial variation of the 
phase profile. 
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Supplementary online material 

Here we provide a theoretical description of the 
model used to describe collective excitations in an 
ensemble of identical meta-molecules. 



Appendix A: Theoretical Model 

We consider an ensemble of N meta-molecules 
placed at positions fj ( j = 1, . . . , N ) arranged in 
a two-dimensional rectangular lattice. This lattice 
is driven by an external, phase modulated, beam 
with electric field Ei n (f, t) as described in the text. 
In general, the charge and current densities for 
the single-molecule oscillatory modes arc described 
by a polarization density Pj(f,t) and magnetiza- 
tion density Mj(f,t). For simplicity, we assume 
that the spatial extent of a meta-molecule is much 
smaller than the wavelength of light and that the 
electric quadrapole moment of the meta-molecule 
vanishes. Each meta-molecule is described by a 
single mode of current oscillation whose contri- 
bution to the scattering of electromagnetic wave 
we then approximate by its electric and magnetic 
dipole moments. The state of the meta-molecule j 
is represented by the dynamic variable Qjit) with 
units of charge and whose spatial characteristics 
have been separated out into time independent 
functions ip(r) and w{f). Explicitly, 

Pi = fj) * dj5(f- fj) (Al) 

Mj = I } (t)w(f- fj) ~ ft&if- fj) (A2) 

where Ij(t) — dQj/dt is the current. In this work, 
we assume the electric dipoles are aligned with 
the incident electric field and that the magnetic 
dipoles are oriented along the driving field's prop- 
agation direction so that only the electric dipoles 
are driven by the incident beam. 

The coupled dynamic equations for the current 
oscillations and electromagnetic fields can be sys- 
tematically derived from the Lagrangian formal- 
ism and will be presented elsewhere [14j. In the 
absence of radiation losses, each current oscillation 
behaves as a simple LC circuit with self capaci- 
tance C and self inductance L. The emitted radi- 
ation, results in a damping of the current oscilla- 
tion, which we treat here using the phenomenolog- 
ical line- widths Te and Tm associated with elec- 
tric and magnetic dipole radiation respectively. In 
the multi-molecule system, however, the radiated 
fields couple to the dipole moments of all other 
molecules in the system. It is convenient to rep- 
resent the current oscillations in terms the dimen- 
sionless normal variables 

w-5ss(% +i %)- (A3) 



where $j = J d 3 rB ■ w describes the effective mag- 
netic flux through the meta-molecule j, and C and 
L are the effective self capacitance and self induc- 
tance of a single meta-molecule. When the evo- 
lution of the input driving fields and the meta- 
molecules themselves have frequency components 
in a narrow bandwidth Sfl about the resonance 
frequency loq = \j\fLC such that 8^1 /ujq <C 1, 
we may neglect any frequency dependence spon- 
taneous emission rates or coupling coefficients be- 
tween the molecules. Further, assuming that r^, 
Tm, 6fl <C Woj we may make the rotating wave 
approximation, in which we neglect terms oscillat- 
ing at high frequencies ~ 2ll>q. In these limits, the 
equations of motion are 

b(t)=Cb(t) + F(t), (A4) 

where b = (b\, . . . , &at) t , the coupling matrix 

C = - Te +™ 1 + ^\ (TeGe + T M g M ) 

+ \y/T E T M {G X +GZ), (A5) 
the driving field contribution is given by 

e^ ot 8 in (t) , , , 

F = i-=^±J-, (A6) 
V2w y/L 

and £i n is the electromotive force induced by the 
driving field £j = J d 3 r E in (f, Q) ■ ip(r — fj), The 
dimensionless coupling matrices Ge, and Gm re- 
sult from interactions with the electric or mag- 
netic fields scattered from the electric or magnetic 
dipoles respectively, while the matrix Gx accounts 
for the electric (magnetic) fields produced by the 
magnetic (electric) dipoles. Equation (|A4[) corre- 
sponds to the integral representation of Maxwell's 
wave equations, as described by Eq.(l) in the 
text, and can be efficiently solved as a linear sys- 
tem. In the limits we've considered, there exist 
as many distinct collective eigenmodes of oscilla- 
tion as there are meta-molecules in the system. 
Each collective mode corresponds to an eigenvec- 
tor of the matrix C and has a distinct resonance 
frequency and decay rate, determined by the imag- 
inary and real parts respectively of the correspond- 
ing eigenvalue. In the context of this simplified 
model, the scheme we proposed in this manuscript 
boils down to choosing a spatial profile and detun- 
ing of the driving field so as to excite a linear com- 
bination of these eigenmodes such that the steady 
state solution is localized on the meta-molecular 
crystal. 



Appendix B: Effect of collective modes 

The effect of strong interactions between differ- 
ent meta-molecules and the resulting collective ra- 
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FIG. 4: First row: The four eigenstates of a 2 x 2 
metamaterial system cj>j (j = 1,...,4), respectively. 
The + (— ) sign refers to the positive (negative) ampli- 
tude of the current oscillation mode. The second row: 
a superposition state $ = ^ . | </>.,• dominantly exciting 
one of the dipoles. Here the sites marked by x are not 
excited. 

diative modes can be illustrated by a simple ex- 
ample of a 2 x 2 metamaterial array. If the inter- 



dipolc separation is much smaller than the wave- 
length of light, the resonance frequencies and the 
radiative damping rates are very different from the 
ones of an isolated dipole and they are affected by 
recurrent scattering processes in which the wave 
is scattered several times by the same dipole. The 
four energy eigenmodes of the system can be repre- 
sented as shown in Fig. 2) Here the positive (neg- 
ative) amplitude of the current oscillation mode is 
displayed as + (— ) sign. 

If an incident field can be prepared in such a 
way that it excites an equal superposition of the 
four eigenmodes $ = J2j h^ji then one can im- 
mediately observe the approximate cancelation of 
the optical excitation energy in three of the four 
dipoles in the 2x2 metamaterial array. This sim- 
ple example illustrates the localization in one of 
the dipoles. 



